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ABSTRACT 

AGN jets carry more than sufficient energy to stave off catastrophic cooling of the intraclustcr 
medium (ICM) in the cores of cool-core clusters. However, in order to prevent catastrophic cooling, 
the ICM must be heated in a near-isotropic fashion and narrow bipolar jets are inefficient at heating 
the gas in the transverse direction to the jets. We argue that due to existent conditions in cluster 
cores, the SMBHs will, in addition to accreting gas via radiatively inefficient flows, experience short 
stochastic episodes of enhanced accretion via thin discs. In general, the orientation of these accretion 
discs will be misaligned with the spin axis of the black holes and the ensuing torques will cause the 
black hole's spin axis (and therefore, the jet axis) to slew and rapidly change direction. This model 
not only explains recent observations showing successive generations of jet- lobes-bubbles in individual 
cool-core clusters that are offset from each other in the angular direction with respect to the cluster 
center, but also shows that AGN jets can heat the cluster core nearly isotropically on the gas cooling 
timescale. One implication of our proposed model is that since SMBHs that host thin accretion discs 
will manifest as quasars, we predicts that roughly 1-2 rich clusters within z < 0.5 should have quasars 
at their centers. Also, recurrent accretion via misaligned accretion discs implies that as a population, 
the SMBHs at the centers of cool-core clusters should be spinning slowly. Our model, in fact, requires 
SMBHs to be spinning slowly. Torques from misaligned discs are ineffective at tilting rapidly spinning 
black holes by more a few degrees whereas slowly spinning SMBHs can, under optimal conditions, 
slew by as much as ~ 30° during any one accretion event. 

Subject headings: galaxies: clusters: general cooling flows galaxies: active - X-rays: galaxies: clusters 



1. INTRODUCTION 

The cool-core conundrum poses a critical challenge for 
theoretical models seeking to explain the observed prop- 
erties of clusters of galaxies. The intracluster medium 
(ICM) in cool-core clusters should, according to sim- 
ple theoretical arguments and di rect observations of X- 
ray emission (|Fabian et al.lll984l ). be cooling and drop- 
ping out at prodigious rates, yet only relatively mea- 
ger amounts of cool and cold gas is seen in cluster cen- 
ters. High resolution X-ray and radio observations of the 
cores of these clusters indicate that the dominant central 
galaxy (hereafter referred to as the BCG — the bright- 
est cluster galaxy) invariably shows evidence of active 
galactic nuclei (AGN) behaviour, often in the form of 
powerful bipolar jets and pairs of approximately spher- 
ical depressions in the X-ray emissions, typically inter- 
preted as being due to bubbles of relativistic plasma that 
have been inflated in the ICM by the jets. The inferred 
power of the jets is comparable to the radiative losses 
in the ICM and the case — based on theoretical argu- 
ments and observational inference — in favour of this 
radio-mode AGN energy injection into the ICM being 
the most likel y explanation for the dimini shed cooling, is 
fairly robust (jMcNamara fc Nulsenl 120071) . However, the 
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precise manner in which this energy is injected into and 
distributed within the ICM remains an open question. 
The mechanics of how an apparently narrow bipolar out- 
flow is able to heat the I CM in a near-isotropic fash ion 
is particularly vexing (c.f. iVernaleo fc Reynolds! 120061) . 

Recent detailed studies of individual cool core groups 
and clusters offer a clue to resolving this "isotropy" prob- 
lem. Many of the systems show evidence of multiple 
generations of radio jets and cocoons as well as X-ra;y 



cavity pairs (cflMcNamara et al.J[2001 



Wise et al. Il2007t iForman et al.ll2007 



Dunn et al 



A-ray 
112001 



Giacintucci et al.l 



20111: IRandall et all 120111 : IQSullivan et al. 112012ft . In a 



number of these systems (e.g. M87 [Virgo], NGC 1275 
[Perseus], NGC 4636, NGC 5044, NGC 5813, Hydra A, 
MACS J0913.7+4056 [also known as CL09104+4109 and 
hereafter referred to as CL09] ) , successive generations of 
jet-lobe-bubbles are significantly offset from each other 
in the angular direction on the sky, with respect to the 
cluster center. If the AGN jet byproducts can trace a 
nearly isotropic angular distribution about the cluster 
center, then so should the associated heating. And if 
the timescale over which the jet axis changes direction is 
shorter than the cooling time, the isotropy problem is a 
non-issue. The question then is: What is the most likely 
explanation for the wandering nature of jet axis given 
the conditions in the cores of cool-core clusters? 
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Potential explanations can be classified into two broad 
categories: The first invoke interactions between the 
jet and the "ICM weather", i. e. large-scak{3 turbulence, 
wakes and bulk velo cities induced by mergers and or- 
biting substructure (jHeinz et al.l 120061: IMorsonv et aO 
120101) . The simulation study by IMorsonv et all (|2010D 
illustrates that strong large-scale turbulent flows can re- 
sult in large angular displacement of the low density 
jet material from the jet axis. Undoubtedly, jet-ICM 
interactions must be occurring at some level; however, 
there are a number of factors that suggest that its im- 
pact in cores of typical co ol-core clusters i s less than 
in the simulations. First, IMorsonv et al.l (|2pi0) only 
consider an unrelaxed cluster. However, the X-ray im- 
ages of the core regions of most cool-core clusters do 
not typically show any evidence of disturbances asso- 
ciated with significant ongoing or recent mergers. Ac- 
cording to detailed simula tion studies of the impact of 
mergers on clusters cores (jPoole et al.l 120071 ). a cluster 
with a dynamically unrelaxed core and no other visible 
morphological merger signatures would need to have ex- 
perienced a major merger (i.e. greater than a 1:3 mass 
ratio) within a relatively narrow window approximately 
3.5 — 4 Gyrs ago. It seems somewhat contrived to as- 
sert that all systems with offset jct-lobc-bubblcs fall in 
this category. Moreover, compar ative analyses of X- 
ray and lensing cl uster data {e.g. iMahdavi et al.l [20081 : 
IZhang et al.l 120101) indicate that the cluster cores are, 
by and large, dyn amically relaxed. Add itionally, while 
the simulations of Morson v et al.l ([2010D show that the 
direction of the high velocity jet flow on small scales 
is unaffected by turbulence, the jets in both Virgo and 
Perseus, for example, show evidence of change in orien- 
tation on subkiloparsec to kiloparsec scales (c.f. section 
[H Table 1), as do the jets in a number of Seyfert galax- 
ies (c.f. IGallimore et al.l 120061 and references therein) 
and in CL09 (|OSullivan et al.Tl2012f ). The latter system 
is particularly interesting in that there are clearly two 
very distinct jets: an older jet that stretches linearly out 
to scales of ~ 60 kpc on either side with no evidence 
of bends or distortions of the kind one would expect if 
there was significant turbulence in the ICM, and a sec- 
ond much smaller jet whose axis is offset from the first 
by ~ 17°. 

The alternate class of models invoke changes in the ori- 
entation of the spin axis of the supermassive black holes 
(SMBHs) that are powering the jets. This explanation 
is premised on the understanding that (a) the angular 
position of any particular feature (jet/cavity/bubble) is 
indicative of the direction of the jet axis at the time when 
the feature was formed, and (b) the jet axis is always 
coincident with the black hole's spin axis. The spin ori- 
entation of black holes can change as a result of pre- 
cession (c.f. IPizzolato fc Sokerll2lT05tllLodato fc Pringld 



I21M ILiu fc Chenl 120071: ICamp airelli et al. 1 120074 an j 
references therein) , spin flips (c.f. Merritt fc Ekers Il2002t 
iCampanelli et al.ll2007b1: iKesden et al. 1120101 and refer- 
ences therein) or slewing (or tilting) of the black hole 
spin axis (c.f. iScheuer fc Feilerl [l99a King et aT~l l2005l: 
ILodato fc Pringld I2006D ." We will consider these option 

1 We use the label "large-scale" to distinguish between veloc- 
ity fluctuations induced by mergers as opposed to "small-scale" 
turbulence generated by, for example, supernova blastwaves. 



in more detail in section [51 after we review some of the 
most compelling observations of repeated re-orientation 
of the jet axes in galaxy groups and clusters. For now, we 
will simply note that the only two viable scenarios are: 
(a) precession associated with binary black hole systems 
in the cluster cores, or (b) stochastic slewing of the black 
hole spin axis due to torques from recurring, short-lived, 
misaligned thin discs. In this paper, we argue that the 
observations are best understood in the context of the 
latter model. 

In the next section, we will first review the most 
compelling observational evidence for recurring re- 
orientation of the jet axes in cluster environments. We 
then assess the various scenarios involving changes in the 
direction of the black hole's spin axis and identify those 
that are both plausible and compatible with the obser- 
vations. In Sj3[ we motivate the case for the our preferred 
model, which involves rapid, stochastic reorientations of 
the SMBH's spin axis on relatively short timescales, sum- 
marize the associated physics, and discuss why we expect 
this to be recurrent phenomenon in cores of cool-core 
clusters. In 21 we discuss the astrophysical implications 
of our model. 

2. WANDERING JET AXIS: 
ASSESSING THE OBSERVATIONAL EVIDENCE 

Observations of " Winged" and "X-shaped " radio 
galaxies (XRGs) (c.f. iHodges-Kluck et alJl2010L and ref- 
erences therein) and changes in the orientation of the 
jets on sub-kiloparse c scales in Seyfert galaxies (c.f. 
IGallimore et al.|[2006l and references therein) have long 
hinted at the possibility that AGN jets undergo dis- 
creet changes in direction by large angles on relatively 
short (a few Myrs to several 10s of Myrs) timescales. 
However, recent radio and X-ray studies of AGN activ- 
ity in cool-core groups and clusters offer arguably the 
most compelling evidence for "jet reorientation" . This 
evidence is manifest in a variety of forms ranging from 
large-angle bends in the radio jets, to abrupt changes in 
direction from one jet episode to the next, to multiple 
generations of radio cocoons and X-ray cavities, some 
of which appear to trace out a nearly isotropic distribu- 
tion in projection about the cluster center. Examples 
of such systems include CL09, M87 [Virgo], NGC 1275 
[Perse us], NGC 4636, NGC 5044, NGC 5813, and Hy- 
dra A (lKleinlll999l : lDunn et al.ll200d iForman et al.ll2007t 
Wise et al. ll2007t iGiacintucci et"aLll2011l : iRandall et all 
201lHOSullivan et al. II2012D . 



In Table [TJ we focus on the Virgo and the Perseus clus- 
ters, two systems that have been the subject of a detailed, 
multi- wavelength observational campaign over the course 
of the past decade, and catalogue most of the pertinent 
features seen in the radio and X-ray maps of these two 
systems typically associated with the radio-mode AGN 
activity of the central SMBH. Collectively, these observa- 
tions of misaligned active and relic jets, as well as radio 
lobes, bubbles, and distinct X-ray cavities separated by 
large angles, strongly indicate that the underlying out- 
flow responsible for forming these features changes di- 
rection between energetic jet "episodes" over timescales 

ra nging from few M yrs to few 10s of Myrs. (Fi g. 1 

in iDunn et al] 120061 and Fig. 5 in iForman et ail 120071 
show this visually.) On the sub-parsec and parsec scales, 
this is indicated by changes in the direction of the jets 
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TABLE 1 

Features Indicating Jet Reorientation in Perseus and Virgo 



Feature 



P.A. (North) 3 



P.A. (South) 3 



Timescale 



Fief. 



Comment 



Perseus (NGC 1275, 3C 84) 

Milliarcsecond jet 



170° 



1.2 



Cocoon-like feature at r 
1.2 mas 

Parsec-scale cocoon at r 
12 mas 
Inner bubble 



Outer bubble 
Ghost bubble 



5° - 10° 
~ 345° 

~ 5° 
~ 305° 



' 210° 
' 180° 

- 155° 

- 215° 
< 170° 



1,2,3 

~ 25 yrs 1,3,4 
~ 10 Myrs 5 

15 - 20 Myrs 5 
70 - 75 Myrs 5 



Ancient bubble 
Virgo (M87, 3C 274) 

Jet 



■ 345° 



290° 



- 230° 



115° 



100 Myrs 



6, 7 



Inner lobes (in radio); 
jet/countcr-jct cavities (in 
X-ray) at ~ 2.5kpc 
Series of buoyant bubbles 
in the southeast 

Intermediate East- West 
jets and lobes 



270° - 275° 



270° 



120° - 130° 



6, 7, 8 



160°, 140°, 120°, 100° (At ago ) ~ 6 Myrs 8 

90° ~ 20 Myrs 9, 10 



Outer X-ray cavity 



35° 



■ 40° 



70 Myrs 



IKrichbaum et al.l (1993) 
argue that the jet is nearly 
aligned with line of sight 
(los) very close to the core. 
On kpc scales, however, 
it is ~ 45° relative to 
los, implying a significant 
change in direction. 



unlike here, 
IDunn et al.l 
measured in 



the clockwise direction. 



Based on the shape of the 
northern bubble and the 
kinematics of the trailing 
Ha fi laments, IDunn et al.l 
( 2006) argue that the bub- 
ble must be moving nearly 
perpendicular to the line of 
sight — in the plane of the 
sky. 




Biretta et al. (1995) argue 
that from 0.1 to 1000 
parsecs, the observed jet 
(stretching to the north- 
west) is aligned at ~ 40° 
to los. The counter jet is 
not observed but its direc- 
tion is inferred from the lo- 
cation of the hot spot. 



Ordered from closest to 
furthest in projection 



Following [Klein (1993), we 
identity the Eastern radio 

ar m (Feature C in Fig. 3 

of lOwen et"aH pOOOT )'). ter- 
minating in the "ear-like 
Feature B, and the West- 
ern extension of the jet 
along with "k not D" as jet - 
lobe features. IKleinl {1999) 
refers to Features B and D 
as "intermediate lobes". 

This is th e cavity dis- 
covere d by IForman et al.l 
(2007) and named "outer 
cavity" . 



Outer radio bubbles 



10° 



205° 



100 Myrs 



9, 10 



a P.A. is the angle in the sky at the source location between North and the feature, measured in the counterclockwise direction. 

f llDhawan et a ll (1199811 ■ (2) Lister (2001), (3) Krichb auTri et al.l fL993t) , (4) INagai eTaT] 02003) , (5) IDunn eTTaf! F200&) , (6l lBirettal (1199311 , 
(7 HBiretta et al.l 1199511 . (81 IForman et al.l ( 120071 1. O lOwen et al.l 120001 . ClO lKleinl 1199911 
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themselves, and on kiloparsec scales and beyond, several 
generations of radio bubbles/X-ray cavities are observed 
and their distribution covers almost all projected angles, 
with many cavity pairs orientated nearly collinear with 
the SMBH. These relative changes in the orientation of 
the jet-counterjet features are more consistent with rota- 
tion of the jet axis rather than displacement of the source, 
which wo uld tend to generat e wide-angle tail-like geome- 
tries (e.g. Uetha et"al]|2008l and references therein). 

Conceptually, the most obvious way to affect a change 
in the direction of the BH's spin axis is via accretion. 
Supermassive black holes in BCGs at the centers of cool 
core clusters are thought to be accreting gas from their 
surroundings via geometrically thick, radiatively inef- 
ficient flows at a relatively low rate M < 10~ 3 MEdd 
([Alien et al.ll2006l ). where 



M 



L 



Edd 



Edd 



0.1c 2 



25 M. 9 



M 



vr 



(1) 



is the Eddington mass accretion rate and M.g = 
( M./10 9 M ) is the black hole mass in units of 10 9 M Q . 
In the event that the angular momentum of the geomet- 
rically thick flow is misaligned with the direction of the 
BH's spin axis, the latter will gradually swivel and move 
to align with the former as BH accretes the gas and 
its angular momentum. The corresponding alignment 
timescale is comparable to the black hole's mass-doubling 
timescale (~ M./M). This timescale, however, is far too 
long to explain the observations: For a 10 9 M Q SMBH 
accreting at a rate M sa 0.01 M©yr _1 , the correspond- 
ing timescale is ~ 100 Gyrs. Accretion via geometrically 
thick flow is not a viable option. 

As noted previously, one way to effect a rapid change in 
the direction of an AGN's jet axis is via a spin-flip, which 
refers to the change in the direction of a black hole's spin 
axis as a result of a major BH-BH merger. However, spin- 
flip is an unlikely explanation for the observed features 
in Virgo and Perseus, for instance. One would need to 
invoke several BH-BH mergers over a duration of ^100 
Myrs to account for the 5-7 episodes of changes in the 
direction of the jet axes over this period. This require- 
ment, in and of itself, makes this scenario implausible. 
Additionally, BH-BH mergers typicall y result in progen- 
itors w ith large recoil velocities (e.g., ICampanelli et al71 
l2007a| ) whereas observations suggest that the origin point 
of the jets appears to be unchanging and coincident with 
the clustcr/BCG center. Large recoil velocities can be 
avoided if the spins of all SMBHs involved are low. Al- 
ternatively, large recoil velocities can be avoided if the 
BHs' spins and the binary's orbital angular momentum 
are aligned prior to the merger; however, such align- 
ment s also preclude the possibility of larg e-angle spin 
flips ([Bogdanovic. Reynolds, k Millerl 120071 ). 

The observations can also be understood if the 
spin axis of the black hole is precessing ab out a 
fixed a xis. A number o f auth o rs, including IKleinl 
(fl999h: IPizzolato k Sokerl (12005^1 : IDunn et al.l ([20061 ): 
iFalceta-Goncalves et al. \ ( 20101 ) have invoked precession 



of the black hole spin axis to explain the distribution of 
radio and X-ray observations in Virgo and Perseus. Pre- 
cession can occur if the black hole in question is part 
of a binary black hole system whose spins and the or- 
bital angular momentum are misaligned. The preces- 



sion timescale due to spin-orbit coupling of a binary 
black hole with both black holes of comparable mass is 
([Merritt k Milosavlievic 1120051: IKev k Cornish I [20H 



/ A \ 5/2 
*BH, P rcc ~ 2.4 x 10 7 yrs( — J M. 



-3/2 



(2) 



where A is the binary semi-major axis. The precession 
timescale is certainly compatible with the observations 
and it is not inconceivable that the BCGs in cluster cores 
host binary supermassive black hole systems since the 
BCGs are thought to have been built up via mergers, in- 
cluding those involving giant elliptical galaxies, at z < 1. 
Moreover, unless the binary SMBH arc embedded in a 
massive, gaseous accretion disc, the lifetime of the bi- 
nary will be more than long enough to allow for m ultiple 
precession cycles ([Merritt k Milosavlievic II2005D . How- 
ever, given the conditions in cores of cool-core systems, 
we would expect that if one of the black holes is powering 
jets, the other ought to be too. There appears to be no 
evidence of dual jets in the cores of Perseus and Virgo, 
nor are we aware of any additional evidence suggesting 
the existence of binary supermassive black holes in ei- 
ther Perseus or Virgo (or in most cool core clusters, for 
that matter). We also note that multiple rc-oricntations 
cannot be straightforwardly modeled as a periodic phe- 
nomenon, a chara cteristic that pose s a challenge for the 
precession model ([Dunn "et~all[2006l) . 

Precession can also occur if the black hole is sur- 
rounded by a geometrically-thin accretion disc whose 
angular momentum is misaligned with the spin axis of 
the black hole. A number of authors (e.g. IKleinl 119991 : 
IFalceta-Goncalves et al. I I2010D invoke a misaligned ac- 
cretion disc as an explanation for the changing jet direc- 
tions in Perseus and Virgo. This scenario, however, is 
problematic for a number of reasons: First, as we noted 
above, precession is a periodic phenomenon and the ob- 
served re-orientations arc not. Second, as Dunn et al.l 



(|2006l) and IFalceta-Goncalves et al~l ([20100 find, the 



Perseus observations can only be understood if the black 
hole precession angle is large and the BH undergoes sev- 
eral precession cycles over the course of ^100 Myrs. If 
the precession is disc-induced, this presupposes the exis- 
tence of a large, stable, warped, geometrically thin accre- 
tion disc over the course of the past ^100 Myrs. AGNs 
with geometrically thin accretion discs are, however, ex- 
pected to shine as quasars. The central AGNs in Virgo 
and Perseus do not appear to be quasars and neither arc 
the AGNs in vast majority of cool-core clusters. More 
importantly, as emphasized by lLodato k Pringlel (|2006l ) 
and discussed more fully in section [3l a misaligned ac- 
cretion disc will not only induce precessional torques on 
the SMBH; it will also induce torques that will cause the 
black hole spin axis to slew and move towards alignment 
with the total an gular momentum of the blac k hole+disc 
system (see also INataraian k Pringle II1998D . Once the 
alignment is achieved, precession ceases. Correspond- 
ingly, a BH will typically only undergo at the most one 
precession cycle. 

The possibility that a misaligned accretion disc can 
cause a bl ack hole's spin ax i s to ti l t on timescales of tens 
of Myrs dScheuer k Feilerl 119961 : INataraian k Pringle I 
[19981 lLodato k Pringldl2006l forms the basis of our pro- 
posed model. Specifically, we argue that while AGNs at 
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the center of cool-core clusters typically accrete gas at 
a relatively low rate via radiativcly inefficient, geometri- 
cally thick flows , every so often the mass accretion rate 
will spike and give rise to short-lived, geometrically thin 
accretion discs. In order to account for the observations, 
we require that the angular momentum vector of these re- 
curring discs is generally not aligned with the direction of 
the SMBH's spin. These recurring misaligned thin discs 
are the central to our model because only geometrically 
thin structures can cause the BH to slew rapidly; geo- 
metrically thick flo ws do not appear to have the same ef- 
fect on the SMBH dKmF et al. M2005t [ Fragile et al.H 2007t 
iDexter fc Fragilell201lD . 

3. TITLING SUPERMASSIVE BLACK HOLES 

There are a number of questions to address in the con- 
text of our "misaligned accretion disc" model is: Do the 
existent conditions at the centers of cool-core clusters al- 
low for the recurring formation of geometrically thin ac- 
cretion discs? Are these discs likely to be massive enough 
to cause the SMBH's spin axis to change direction sig- 
nificantly on short timcscales? 

3.1. Recurring, Short-lived High Accretion Events in 
Cool-Core Environments 

Apart from the hot diffuse gas, central regions 
of BCGs in cool-core clusters also appear to be 
threade d by a filament a ry ne twork of cold gas (c.f. M87 
[Virgo ] IForman et all (12001 : CL09 IQSullivan etHI 
pOia and references there in); NGC 4696 [Cen- 
taurusl ICrawford et al.1 (120051) : NGC 1275 [Perseus] 
IConselice et al.1 ([20011) ) that generally appear to cover 
a significant fraction of the 4ir steradian about the 
cluster center. Increasingly detailed kinematical stud- 
ies of the filaments suggest that these likely have 
their origins in number of different physical mecha- 
nisms. Many of these mechanisms will also induce short- 
lived accretion events during which the instantaneous 
mass accretion rate onto the central SMBH can ex- 
ceed O.OlMEdd, the threshold accretion rate above which 
the geometrically thick, radiativcly inefficient flow is ex- 
pected to transition to a geometrically thin, di sc flow 
(jEsin. McClintock fe Nar"avan1 H997t I Jested I2005T ). Here 
we briefly review the three likely mechanism s: bub- 
ble wakes (lHatch et al.ll2006t iPope et al.l|2010D, thermal 
insta b ilities (iPizzolato fc, Sokerl l2005al TMcCourt et al.l 
l2012t ISharma et al.l 12012k and c onvergent small-scale 
t urbulent flows (IHobbs et al.ll201ll ). 

IPope et all (|201C0 present a detailed description of the 
mechanisms by which buoyantly-rising jct-inflatcd bub- 
bles transport cool gas out of the BCG cores as they 
detach and start to rise, along with the properties of 
the "cool wake" . For quasi-spherical bubbles of the type 
typically seen in the X-ray observations, this wake is ex- 
pected to involve ~ 10 8 M© of cool gas per bubble. The 
trails of cool gas behind the rising bubbles can stretch 
for several tens of kiloparsecs, with filamentary material 
close to the bubble traveling with the bubble, and away 
from the cluster centre, and the material further down 
the filament (and closer to the cluster centre) falling back 
towards the cluster centre. Given this kinematic pro- 
file, it is expected that every so often fragments of cool 
gas (M ~ 10 5 -10 7 M©) will break away from the fil- 
amentary wakes and drip onto the BH. Since the fila- 



ments are primarily radial in orientation, we expect that 
the fragments will have intrinsically low angular momen- 
tum relative to the SMBH and hence, ultimately settle 
into a subparsec-scale accretion disc. The angular mo- 
mentum of the wake (and therefore, the alignment of 
the accretion disc that forms from material breaking off 
from this wake) will be more or less perpendicular to 
the black hole's spin/jet axis because the bubbles are 
inflated along the axis of the jet that created them and 
buoyantly rise radially outward. Anticipating the discus- 
sion below, each jet/bubble outburst in effect sets into 
motion a chain of events that result in the reorientation 
of the BH's spin axis and hence, the direction of the jet. 

In the case of turbulent gas, iHobbs et al.l (|201lD fo- 
cus on the formation and evolution of dense, low-angular 
momentum, filaments created by convergent turbulent 
flows. In their simulations, they find that such filaments 
involve only a small fraction of the total available gas. 
They find that these streams fall "ballistically" onto the 
SMBH, where the gas is presumed to settle into a sub- 
parsec scale accretion discs. Translating the angular mo- 
mentum of individual filaments crossing the lower bound- 
ary of their simulation s into an or i entatio n vector for 
the accretion disc, the IHobbs et al.l ([20111 ) results sug- 
gest that not only are successive disc orient ations uncor- 
recte d, they can be offset by large angles. IHobbs et al.l 
(|2011l) do not compute the mass spectrum of these accre- 
tion events nor do they estimate the frequency. We note 
that their simulations do not include self-gravity, cooling, 
star formation and stellar feedback, which may explain 
why the filaments do not fragment into clouds. Typically, 
one would expe ct turbulence to give rise to clouds with 
mass spectrum (|Hopkins. Quataert fc Murravll2012l ) 



diVciocMT 1 ' 8 dM ch 



Mci < / g 2 as M gas 



(< R), 



(3) 



where / gas is the fraction of gas mass relative to the total 
enclosed mass. The typical mass of cold gas in cool-core 
cluster BCGs is ~ a few xlO 10 M^ within the central 
lOkpc (lEdgdl2001t [Salome fc Combed 12003 lEdee et a.l 
I2010T ) . which corresponds to / gas ~ 0.1. In this case, the 
number of clouds with mass M c \ > 10 6 M© is ~2,500. 
(We choose this mass threshold because less massive 
clouds have negligible impact on the SMBH; we discuss 
this further in §3.2 below.) Only a fraction of these 
clouds, with specific angular momentum small enough 
such that the circularization radius is ~ 0.1 pc, will give 
rise to accretion discs. Assuming that the velocity distri- 
bution of the clouds is isotropic, with velocity dispersion 
cr c i « 300 km/s, and that the radial distribution of the 
clouds is approximately isothermal (n c \ oc r~ 2 ), the mean 
time between accretion events is ~ a few x 10 yrs. 

Recent nu merical simulation studies show that thermal 
instabilities (ISharma et al.l 120101: IMcCourt etHI 120121: 
ISharma et al. 120121 ). as well as small-scale supersonic 
turbul ence — induced, fo r example, by supernova blast- 
waves (jHobbs et al.ll201lT ) associated with s tar formation 
in th e central regions of the BCGs {e.g. IBildfell et all 
I2008D — can also give rise to gaseous streams and 
high density clouds. The fraction of these structures 
with low angular momentum will rain down upon the 
AGN giving rise to stochastic events of enhanced ac- 
cretion that c an approach the Eddington rate for short 
periods (e.g., ISharma et al.l 120121 : iGaspari et al.l 120121 : 
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iHobbs et al.l 1201 ID . Within the therm al equilibrium 
model for cool cores (|Sharma et al.|[20T2T ). the occurrence 
of cold filaments and blo bs depends on the fee dback heat- 
ing efficiency (Fig. 11 in lSharma et al.ll2012| ). If the effi- 
ciency is high, the core reaches equilibrium at low densi- 
ties and will not be prone to forming thermal instabilities. 
At low efficiencies, however, cold gas is present most of 
the time. In general, we expect the gas in the cores of cool 
core clusters to undergo distinct cooling "events" sepa- 
rated by several core cooling times (i.e. a few 100 Myrs). 
However, an individual cooling "event" will spawn mul- 
tiple cold blobs with slightly different angular momenta, 
which will give rise to closely spaced accretion events, 
and therefore BH realignments, with time lags ~ 10 6 - 
10 7 years, i.e. , of order the free-fall time. These can ex- 
plain angular offsets between bubbles/cavities separated 
by < 100 Myr in age. 

3.2. Misaligned accretion discs and spinning black holes 

Given an accretion disc whose initial rotational axis is 
misaligned with the spin axis of the black hole, frame 
dragging by the rotating black hole will induce a torque 
on the inner regions of the disc that will cause it to pre- 
cess diffe rentially, an effect known as L ense-Thirring pre- 
cession. iBardeen fc Petterson I (|1975l) showed that the 
viscous forces in the disc will damp the differential pre- 
cession, and force the angular momenta of the disc to 
aligi fl with the to t al angular momentum o f the system 
([King et al. I [20051 Hodato fc Pringlel[200l . The align- 
ment of the disc proceeds from inside out, with the tran- 
sition radius between the inner disc and the outer disc 
demarcated by a warp. Since the process is primarily 
driven by frame-dragging, the influence of which falls off 
rapidly with distance, the warp will not propagate out 
across the entire disc. Instead it will stall at a radius 
i? w , where the rate at which the disc is twisted by Lense- 
Thirring precession is balanced by the rate at which vis- 
cous torques can dissipate the twist dScheuer fc Feilerl 
[1991 IKing et al. I [20051 : lLodato fc Pringlel 12006ft The 
warp radius is given by (see Eq. 22 in lKing et al. 1120051 ) 



2.6xl0- 3 j 5/ 1 8 



M, 



-1/4 



0.04 



- 1 / 2 / OLijoty 

V 30 



-1/8 



pc, 



(4) 



where M . 4 = (M/0.04M E dd) and jo.i = (j/0.1). Here, 
j is the black hole spin parameter scaled to J m ^ nax = 

(G M. 2 /c), the maximum angular momentum of a Kerr 
black hole: < j < 1. King and colla borators have ar- 
gued in a series of related articles (c./. IKing et al"1l2008l 
and references therein) that if BHs grow primarily via re- 
current randomly orientated accretion events, they will 
tend to have low spins. For this reason, we have cho- 
sen to scale j to a fiducial value of 0.1. Also. a\ in 
the above relationship is the iShakura fc Sunvaevl ([1973D 
viscosity parameter characterizing the radially outward 
transport of gas angular momentum and inward trans- 
port of matter flow, and cx2 is the viscosity associated 

2 Here, we use the term "align" loosely to refer to either co- 
alignment or counter-alignment. Details of how the disc orientation 
evolves is immaterial to the present discussion. What is important 
is that the gas is subject to torques that drives a change in its 
orientation. 



with vertical motions in the disc. For the purposes at 
hand, we assume that over the regions of interest, the 
ct\ parameter (the usual accretion disc a parameter) is 
approximately constant and adopt u\ ~ 0.1 as a fidu- 
cial value. This is consistent with results from MHD 
simulations of magnetized accretion flows, which find 
that the effec tive value of a± is ~ 0. 1 over the bulk of 
the flow (e.g..lHawlev fc Krolik |[200ll l200l). Moreover, 
lLodato fc Pried ([20101) show that for a strongly warped 
disc, a2 ~ 3. For discs with ot\ > H/R, which should 
hold true for thin disc accretion, the warp will propagate 
through the disc diffusively on timescale a^j o.\ shorter 
than the mass accretion timescale, which is determined 
solely by ol\. Finally, we note that we have scaled the 
mass accretion rate to 0.04MEdd, which corresponds to 
1 Mgyr" 1 for a 10 9 M . 

The disc is not the only structure to experience a 
torque. The black hole too will experience equal and 
opposite torques exerted by each differentially process- 
ing disc annulus, which in turn cause the black hole's 
spin axis to change direction. In the case of a disc whose 
outer radius is less than i? w , the torques on the black 
hole will exist for only a very short time while the warp 
expands out across the disc; such discs do not impact 
the black hole in any significant fashion. For this reason, 
we will restrict ourselves to discs whose radius extends 
beyond R rm w In such cases, the torques on the black 
hole is primarily d ue to the gas flowing through the warp 
(jKing et al.ll2008l ). These torques can be resolved into 
two independent components: one that drives the preces- 
sion of the hole's spin about the total angular momentum 
axis of the disc+BH system and another that drives the 
alignment of the black hole's sp i n wit h the total angular 
momentum. lLodato fc Pringlel (|2006h and iMartin et all 
(|2007| ) have investigated these two processes in detail and 
find that the BH precession timescale is comparable to, 
and in realistic cases perhaps even a factor of a few longer 
([Martin et al.l [20071 ) than, the timescale over which the 
black hole's spin will align with the total angular mo- 



mentum of the BH+disc system: 



^BH 



prcc 



> t 



align i 



where the black hole a lignment timescale (Eq. 15 in 
Lodato fc Pringlel 120061: see also iScheuer fc Feilerl 119961: 



Nataraian fc Pringle Ill99cl : 



Malign « 2.0 x 10 6 i U / 16 



M" 7/8 M ~ 1/16 



Ct2 \ 1 / 4 ( 012 /a 



/02 

V 3 



30 



-15/16 



yrs. (5) 



Consequently, the black hole is unlikely to execute more 
than one precession cycle, if that, and once alignment 
is achieved, the torques driving the precession vanish as 
well. 

Whether the BH attains full alignment during a given 
accretion event depends on the size (mass) of the ac- 
cretion disc. The disc must be sufficiently long-lived to 
sustain gas flow (and therefore, the torques) over the du- 
ration taiign- In other words, 



M dialign > 2.0 x 10 6 Jo 11 / 16 M.9 



15/16 



(¥)' 



/ a / I \ -15/16 

/ 4 / a 2 /ai \ 
\ 30 



M . (6) 



In general, the above scenario for the slewing of the 
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BH's spin axis docs not require the BH to achieve full 
alignment with the total angular momentum during each 
and every accretion event, only that the BH's spin axis 
change direction. However, as we discuss below, tilts in- 
duced by discs with Md < 10 6 Mq are too small to be 
of interest. 

At the same time, the amount of gas that a black hole 
can accrete in the thin-disc mode during any one accre- 
tion episode cannot be arbitrarily large regardless of the 
amount of gas that is channeled into the nuclear region 
during the event. This upper bound arises because accre- 
tion discs are susceptible to gravitational fragmentation 
beyond a radius R{ las where the disc mass exceeds (c.f. 
Goodmanl [200l [Thompson et al.l 120051 King fc Pringli 
20071: King et al.ll2008l and references therein) 



M d (< i? frag ) = / 



M„ 



(7) 



where the factor, / ~ a few, represents the uncer- 
tainty in this relationship due poorly understood de- 
tails such as the extent to which the magnetic pressure, 
radiation pressure and stellar feedback augments ther- 
mal pressure due to viscous heating and contributes to 
the stability of the disc. Adopt ing a simple model of 
iCollin-Souffrin fc Dumont! ()1990l ) to describe the prop- 
erties of the sub-parsec-scale thin accretion disc (see also 
King et alJl2008f Pl. the total disc mass inside radius R. 
where R denotes the distance from the BH, is 



-4/5 



M d «3.1xl0 6 M^ 5 4 M4 /5 (|i 
and the accretion disc scale height H is 



H 



i-3 lOrVS 



— w 1.6 X 10- 5 M c 



3/20 / Oil 



M.q — 

1 0.1 



-1/10 



R 



1/20 
0.05 ' 



(8) 



(9) 



Here, i?o.05 = (i?/0.05pc). Using these equations, we 
find that 

i?ftag « 0.1 f B M < 4 M. g [—) pc. (10) 



Here, f$ = (//5). Most of the gas at R > i?f rag is ex- 
pected to either turn into stars or be expelled by stars 
that do form on timescales much shorter than those that 
govern the accretion flow. Since only gas that flows 
through the war p induces a torque on the black hole 
(jKing et al.l2008D , this means that the maximum amount 
of mass that can participate in the alignment process is 



7.4xl0 6 /f /27 M^ 27 



x M. 23/27 



U.l 



-2/27 



M 



• 



(11) 



Finally, we can estimate the maximum angle i/i max 
through which the black hole's spin axis will slew dur- 
ing any one accretion event where the initial misalign- 
ment between the disc angular momentum and the black 

3 This is a model of a flat steady-state disc and we acknowledge 
that a flat disc and a warped disc are unlikely to have identical 
structures. The use of this model is a necessary simplification in 
order to estimate relevant disc mass and length scales. We note, 
for completeness, that the scaling properties of this d isc model are 
similar to that derived by Shakura & Sunyacv (1973). 



hol e spin vector is c haracterized by angle 6i. Follow- 
ing [King^Fal] (|2008t) . we identity Jd, the "disc angular 
momentum" , as the total angular momentum of the gas 
flowing through the warp. Consequently, if AM gas is the 
total gas mass that flows through the warp during an 
accretion event, then 



J d = (G M. R w ) 1/2 (AA/ gas ) , 



(12) 

where AM gas < M djmax . Straightfoward geometry then 
gives 



Smip max = — Shi(0i - Vmax). 
•Jh 



(13) 



The above equat ion differs slightly from that given by 
King et all (|2008l ) because the latter's derivation implic- 
itly assumes that V'max 1. The extent of the tilt in- 
duced by an accretion disc depends on the initial mis- 
alignment between the disc's angular momentum and the 
BH's spin axis as well as on the ratio 



Jd = V2 / AM,, 



Rw 
Rs 



1/2 



(14) 



where R s = 2GM/c 2 10~ 4 M. 9 pc is the 
Schwarzschild radius. To effect a tilt of ~ 5° in the spin 
axis of a slowly spinning {i.e. j ~ 0.1) SMBH, the min- 
imum amount of gas required is AM gas ~ 10 6 M Q . A 
maximal disc, i.e. AM gas ~ Md jmax , can tilt a a slowly 
rotating SMBH by as much as ~ 30°. For completeness, 
we note that rapidly spinning black holes {i.e. j ~ 0.9) 
are much harder to tilt. A maximal disc can tilt a rapidly 
spinning {j ~ 0.9) supermassive black hole by, at the 
most, — 7°. 

4. DISCUSSION, IMPLICATIONS AND SUMMARY 

To recap, supermassive black holes at the centers of 
cool-core clusters will, in addition of the accretion of 
hot diffuse gas from its surroundings, experience short- 
lived, recurring episodes during which the instantaneous 
mass accretion rate onto the central SMBH can approach 
~ 1 M© yr _1 . We expect this phenomena to be ubiq- 
uitous; it is a byproduct of a number of very different 
and distinct physical processes expected to be operat- 
ing in cores of cool-core clusters. During such episodes, 
the flow in the immediate vicinity of the black hole will 
transition to a geometrically thin flow and establish an 
accretion disc. 

Current generation of simulation studies looking at the 
wide variety of processes at play in cool cluster cores are 
not detailed enough to provide a full statistical charac- 
terization of the spectrum of gas mass involved in indi- 
vidual high density accretion events, the distribution of 
the time between accretion spikes, or the direction of the 
rotation vector of the resultant accretion disc relative to 
the spin axis of the SMBH. Nonetheless, the simulations 
do suggest that all of these will vary stochastically from 
event to event and specifically, the orientations of suc- 
cessive discs will be uncorrelated. Simulations as well as 
theoretical estimates also suggest that individual accre- 
tion events will involve 10 5 -10 M Q and that the average 
time between the events is ~ 10 Myrs, with the caveat 
that the time between pairs of consecutive events can 
vary widely. This timcscale is in agreement with that 
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implied by the history of AGN activity in Perseus and 
Virgo, as summarized in Table 1: Perseus and Virgo data 
both suggest ~ 5 events over the past 10 s years. 

When an accretion event gives rise to an accretion disc 
whose orientation is misaligned with the spin axis of the 
black hole, the resulting torques between the gas flow- 
ing through the disc and the black hole will cause the 
latter's axis to swivel and change direction. During any 
one event, the magnitude of the change in the direction 
of the BH's spin axis depends on the degree of the ini- 
tial misalignment between the disc's angular momentum 
and the BH's spin axis, the magnitude of the disc angu- 
lar momentum, and the magnitude of the SMBH's spin 
angular momentum. We find that under optimal condi- 
tions, a slowly rotating SMBH's spin axis can slew by as 
much as ~ 30° . In general, the direction of the spin axis 
after an accretion event will be uncorrelatcd with the di- 
rection of the old spin axis and over the course of many 
such events, the orientation of the black hole's spin axis 
will execute a random-walk over Air steradian. 

Drawing on the considerable body of work indicat- 
ing that jet production is most efficient when the 
accretion f low is geom e trical l y thick and suppressed 
otherwise (iLivio et al. I 119991 iMeier I 120011 ) (c./.also 
iBenson fc Babull 120091 and references therein), we ex- 
pect that the jets will briefly wane in power during the 
thin disc phascQ But once the accretion disc drains away 
and the geometrically thick flow re-establishes, the jets 
will also resume. Since the jet axis is the same as the spin 
axis, the post-tilt jets will point in a different direction 
from the pre-tilt jets. This scenario offers the simplest ex- 
planation of, for example, the apparently abrupt change 
in direction between an older (larger) relic jet in CL 09 
and a subsequent jet episode (jOSullivan et al. II2012I ). 

The mechanism/model outlined in this paper and sum- 
marized above has a number of important astrophysical 
implications: 

• The recurring reorientation of AGN jets due to tilt- 
ing of the SMBH spin axis will naturally give rise to 
distinct, randomly oriented jets/lobes/cavities ob- 
served in cool-core clusters such as Perseus, Virgo 
and CL09. 

• If the jets typically change directions several times 
over the cooling time in cores of cool-core clusters 

— as we expect them to, since the latter is typically 
a ~ 100 Myrs — the resultant heating will, like 
the cavities and bubbles, also be distributed over 
a large portion of Air steradians. We suggest that 
this process plays a significant role in the resolution 
of the "isotropic heating" puzzle. 

• Since geometrically thin accretion discs are radia- 
tivcly efficient, our model predicts that whenever 
one forms, the host AGN ought to transform into a 
quasar. The quasar will be short-lived because the 
discs do not involve a lot of mass and are expected 
to drain quickly. We can estimate the likelihood 

4 This statement should not be interpreted as our suggesting that 
jet activity cannot arise during the quasar phase. We appreciate 
that at sufficiently high mass accretion rate, the flow should be 
able to sustain both behaviours. However, we do not expect such 
configurations to be common in cluster environments. 



of catching an AGN "in the act" as follows: The 
Perseus and Virgo data both suggest ~ 5 events 
over the past 10 8 yrs. Assuming that a typical ac- 
cretion event involves ~ 10 6 M Q of gas and that 
the lifetime of the resulting disc is ~ Myrs, we ex- 
pect a quasar duty cycle of ~ 5%. The rest of 
the time, the AGNs will be observed in the radio- 
mode, driving jets. The number of rich clusters 
(T x > 2 keV) within z = 0.5 that have X-ray ob- 
servation s is ~ 250. Approximately 35% of rich 
clusters (jEckert. Molendi fe PaltaniH2011l ) tend to 
be strong cool core systems and of these, 35% show 
evidence of sig nificant multiphase gas component 
in their cores ((McDonald et al.l 120101) . Combin- 
ing these, we would expect 1-2 of the z < 0.5 
rich clusters to be in the quasar phase. In fact, 
there are two: MA CS J0913. 7+4056 [also known 
as CL09104+4109] (jOSullivan et al. ll20ll at z = 
0.44 which hosts a dust enshrouded type 2 QSO 
at its centre, and a z = 0.3 cluster that hosts 
a highly luminous ra dio-quiet quasar, H1821+643 
(|Russell et al. If2010h . Both quasars are located at 
the centres of cool-core clusters. 

• The manner in which SMBHs acquire most 
of their mass has an impor tant bearing on 
their spin distribution (c.f . iBerti fc: VolonterU 
l200l iFanidakis etHI 12011 and in turn, on 
their jet efficiencies (c.f. iNemmen et al. I [20071 
IBenson fc Babull [20 09). If BHs grow primarily via 
recurrent accretion events of random angular mo- 
menta, then as King and collaborato rs have ar- 
gued in a series of related articles (c.f. King et all 
I2008L and references therein), they will tend to 
have low spins. The alternative model assumes 
the b lack holes grow via c oherent accretion flows 
(e.g. IVolonteri et al.l [2005D and predicts that the 
bulk of the SMBHs are rapidly rotating, with 
j « 0.9 whe n the effects of the je ts are taken 
into account (iB enson fc Babull 120091 ). However, 
even in such models, once the host dark matter 
halo of the SMBH grows to cluster-scale, i.e. 10 14 - 
10 15 M Q , and the system itself evolves into a cool 
core cluster of the kind discussed in this paper, 
the subsequent evolution of the SMBH ought to be 
dominated by d i screet stochastic accretion events. 
iVikhlinin et al.l (|2007l ). having found a sharp de- 
cline in the number of cool core clusters at z > 0.5 
within sample of X-ray selected clusters, suggest 
that cool-core clusters are a relatively recent phe- 
nomenon. Assuming that this is true and that 
these systems formed about 5 Gyrs ago, that the 
mean accretion rate over this period is comparable 
to that suggested by the Perseus and Virgo data, 
and that typical accretion events involve ~ 10 6 M Q 
of gas, we find that the total gas mass accreted via 
the stochastic process is AM gaSjacc ~ 2.5 x 10 8 M , 
which is a substantial fraction of a SMBH's mass. 
This, in turn, suggests that the SMBHs at the cen- 
tres of cool-core clusters ought to be slow rotators. 
If the cool core clusters arc older than 5 Gyrs, 
as suggested by r ecent SZ detection s of cool-core 
clusters at high z dSemler et al.ll2012T ). the case for 
slowly spinning SMBHs is even stronger. 



Rapidly Reorienting AGN Jets in Galaxy Cluster Cores 



9 



• Finally, we note that slowly rotating SMBHs arc 
not only the result of the model described in this 
paper, our model requires SMBHs at the centers of 
cool-core clusters to be slowly rotating in order to 
explain the observed angular offsets between suc- 
cessive generations of jet-lobe-bubbles. If SMBHs 
in such environments are shown to be spinning 
rapidly, our proposed mechanism cannot explain 
these observations. Thin discs of the kind described 
in this paper cannot tilt rapidly spinning SMBHs 
by more than a few degrees. 
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